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Silicon-based quantum cascade lasers (QCLs) offer the prospect of integrating coherent THz ra- 
diation sources with silicon microelectronics. Theoretical studies have proposed a variety of n-type 
SiGe-based heterostructures as design candidates, however the optimal material configuration re- 
mains unclear. In this work, an optimization algorithm is used to design equivalent THz QCLs 
in three recently-proposed configurations [(001) Ge/GeSi, (001) Si/SiGe and (111) Si/SiGe], with 
emission frequencies of 3 and 4 THz. A systematic comparison of the electronic and optical proper- 
ties is presented. A semi-classical electron transport simulation is used to model the charge carrier 
dynamics and calculate the peak gain, the corresponding current density and the maximum oper- 
ating temperature. It is shown that (001) Ge/GeSi structures yield the best simulated performance 
at both emission frequencies. 
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I. INTRODUCTION 

Terahertz quantum cascade lasers (THz QCLs) are 
semiconductor devices in which electrons are transported 
through a periodic multiple quantum well heterostruc- 
ture, with a radiative transition in each period.— THz 
QCLs have numerous potential applications including ra- 
diation sources for medical and security imaging, and lo- 
cal oscillators in astronomy and remote gas sensing^— 

All THz QCLs to date have been fabricated from III 
V compound semiconductors (for example Ref. and 
8). However, Si-based QCLs could offer a number of 
significant advantages. Mature Si processing technol- 
ogy may reduce costs and allow integration with con- 
ventional electronic devices. Existing THz QCLs op- 
erate only at cryogenic temperatures (currently below 
186 K for resonant-phonon QCLsl or 116 K for bound- 
to-continuumS) , but the high thermal conductivity of Si- 
based structures could enable heat to be extracted more 
effectively, and hence allow higher operating tempera- 
tures. III-V QCLs are also limited to THz emission at 
frequencies lower than 4.9 THz^S owing to the strong ab- 
sorption in the Reststrahlen band. However, this limita- 
tion does not exist in non-polar group-IV materials. 

Although the indirect bandgap in Si has, so far, frus- 
trated efforts to develop an interband laser, this is not 
an issue for intersubband devices such as QCLs because 
the radiative transitions occur between subbands within 
the same valley of an energy band. Mid-infraredii and 
THz 12 intersubband electroluminescence has been ob- 
served fromp-type SiGe/Si quantum cascade structures. 
However, dispersion relations for holes in these struc- 
tures are quite complicated, owing to the contributions 
from multiple valence bands, and in recent years, atten- 
tion has switched toward n-type structures. This greatly 
simplifies the device design process and may result in a 
lower spectral linewidth (and hence greater peak gain) 
than that of p-type structures. Early design proposals 



included n-type Si/SiGe structures that exploited transi- 
tions in the A valleys of (001) Si quantum wells (QWs), 
with SiGe barriers. 13 However, a range of alternative ma- 
terial configurations have been considered in recent years. 
L valley transitions in (001) oriented Ge/GeSi have at- 
tracted greatest attention,— ~— although transitions in 
the A valleys of (111) oriented Si/SiGe^^ the T val- 
ley of Ge/GeSi^ and the L valleys of Ge/GeSiSr3 have 
also been considered. 



Many properties relating to the bandstructure and car- 
rier transport have a strong effect upon the gain of QCLs. 
To date, no quantitative comparison has been made of 
the expected performance of THz QCLs in different Si- 
based material configurations. In this paper, we present 
such a comparison by simulating equivalent devices in 
the (001) Si/SiGe, (111) Si/SiGe and (001) Ge/GeSi ma- 
terial configurations, using a detailed semi-classical rate- 
equation approach. Devices emitting near 3 and 4 THz 
were designed for each material configuration by rescal- 
ing a recent bound-to-continuum design 14 - according to 
the effective mass of the material and then applying an 
automated design-optimization algorithm. In section ITTl 
we calculate the conduction band potentials for each ma- 
terial configuration and the range of energies within QWs 
that can be used for QCL design. Section Mil describes 
our model for simulating carrier transport, gain, and cur- 
rent density in QCLs. The design-optimization algorithm 
is described in section llVl and a summary of the optimized 
devices is presented. Finally, section fVl presents a com- 
parison of the simulated performance of devices in each 
material system. 
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TABLE I: Material parameters for Si and Ge. 
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constant of a thin Sii-^Ge^ layer was found using 
a(x) = a Si (l - x) + a Gc x - b how x(l - x), 



(4) 



where &bow = 0.2733 pm2i is a bowing factor^ The layer 
was assumed to deform uniformly, such that the lattice 
constant matched that of a thick substrate material, a s . 
The resulting strain in the plane of epitaxial growth was 
given by e\\ = (a s - a) /a. 

The hydrostatic deformation was found for each set of 
conduction band valleys using 
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(5) 



where E<i+ ^u~- ctv is the bandgap deformation potential 
for the A or I valleys and Tr e' is the trace of the strain 
tensor, where^i 



II. CONDUCTION BAND STATES 

A. Model solid approximation 

The model solid approximation— was used to calculate 
the conduction band offset between a strained Sii_ x Gea; 
alloy and a substrate material. The difference between 
the average of the light-hole, heavy-hole, and spin-orbit 
split off valence band edges in the two materials was used 
as a reference energy, as it is almost independent of strain 
and crystal orientation,— The value of this property was 
interpolated from empirical pseudopotential data as 26 



AE V = (0.47 - 0.06a; s )(:z - x s ), 



(1) 



where x s is the Ge fraction in the substrate. The valence 
band maximum in an unstrained bulk alloy was found 
using 



E, 



bulk 



AE„ 



1 



-A 



so, 



(2) 



where Aso is the spin-orbit splitting energy. Material pa- 
rameters for alloys were found by interpolating from the 
bulk Si and Ge values in table HI These parameters yield 
an offset of 0.55 eV between the top of the valence bands 
in bulk Ge and Si, which matches a recently-measured 
value for weakly-strained Ge and Si filmsiSS Very similar 
parameters have also been shown to yield close agree- 
ment with experimental measurements of intersubband 
absorption energies in Ge/GeSi QWs<2& 

The low-temperature indirect bandgaps for the A and 
L valleys in an unstrained alloy (in eV) were taken as 2 - 



E* = 1.155 - 0.43a: + 0.0206^ 
E% = 2.010- 1.270a;, 



(3) 



and the unstrained conduction band edge for a given val- 
ley was found using E^ ulk = E, 



bulk 



E n 



The effects of hydrostatic and uniaxial strain on the 
band edge were determined as follows. First, the lattice 
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for (001) and (111) oriented epilayers, respectively. In the 
above equations, Cn,Ci2, and C44 are elastic constants. 

Uniaxial strain leads to splitting of the A valley de- 
generacy in (001) oriented layers. The energy shifts are 
given by^i 
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for the valleys with their major axes perpendicular and 
parallel to the growth-direction respectively, where 
is the uniaxial deformation potential for the A valleys. 
Similarly, L valleys in (111) oriented layers are shifted by 
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for the valley with its major axis in the growth direction, 
and the three other valleys respectively. Uniaxial strain 
has no effect upon the A valleys in (111) layers, or the L 
valleys in (001) layers, owing to symmetry. 

Finally, the energy of a given conduction band mini- 
mum in a strained layer was found relative to the average 
substrate valence band using 
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B. Available energy range 



In this section, we determine the energy ranges within 
jroup-IV QWs that can be reliably exploited for THz 



1.1 



1.05 



> 



1- 



&0.95 

U 

w 0.9 
0.85h 
0.8, 



Available 
energy range 



A , 

— L 







4 6 
Growth axis [nm] 



FIG. 



1: (Color online) Spatial variation of A and L 
conduction band minima in a 
Sio.15Geo.85/Ge/Sio.15Geo. 85 heterostructure on a 
Sio.025Geo.975 virtual substrate. The energy range 
available for QCL design is defined here as the region 
between the bottom of the QW and the next-lowest 
conduction band valley. 



QCL design. It is insufficient to simply calculate the 
depth of a QW as there are multiple conduction band 
valleys within the energy range of interest. This can 
degrade device performance by introducing undesirable 
intervalley scattering processes. To avoid this problem, 
we consider the energy difference between the bottom 
of the well, and the next-lowest conduction band min- 
imum (which may be in either the well or the bar- 
rier). This is illustrated in Fig. [1] for the case of a 
Sio.15Geo.s5/Ge/Sio.15Geo. 85 QW on a Sio.025Geo.975 sub- 
strate. Here, the L valleys form the bottom of the well 
and the usable energy range is limited by the A2 valley 
minima in the barriers. 

Fig. [5] shows the conduction band minima in QWs that 
consist of either a Si or Ge well surrounded by SiGe bar- 
riers. The energies of the valley minima in the well and 
barriers are plotted as a function of the barrier alloy com- 
position. In each case, the energies are expressed relative 
to the bottom of the QW, and usable energy range is 
therefore given by the lowest line in the plot. 

Fig. [ ^a)| shows the results for a (OOl)-oriented 
SiGe/Si/SiGe QW. Here, the A 2 valleys form the bot- 
tom of the well in the Si layer. The total depth of the 
QW is given by the energy difference between the A2 
minima in the barrier and the well. As we shall see in 
section IIV1 the barriers in (001)-Si/SiGe QCL designs 
may need to be thinner than 1 nm, owing to the large 
A2 quantization effective mass. It is, therefore, neces- 
sary to limit the barrier Ge fraction to obtain a lower A2 
confinement potential, and hence a realistically wide bar- 
rier layer. We selected a Sio.sGeo.2 alloy for the barriers, 
which provides a A2 band offset of 95meV. A ^4.5% Ge 
virtual substrate is required for mechanical stability in 



TABLE II: Quantization and two-dimensional 
density-of-states effective masses of conduction band 
valleys in (001) Si, (111) SLand (001) Ge films, 
calculated as described in Ref. |32|. Masses are expressed 
relative to the rest mass of a free electron. 
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our (001)-Si/SiGe QCL designs because the Si wells are 
considerably thicker than the SiGe barrier layers. This 
induces only a relatively low uniaxial strain in the QW 
layers and as a result the A4 minima in the barrier layers 
lie only 5 meV above the bottom of the well. It is, there- 
fore, impossible to avoid the presence of A4 subbands 
within the energy range of interest for (001) Si/SiGc QCL 
designs. 

Fig. Qth)] shows the calculated valley minima for (111)- 
oriented SiGe/Si/SiGe QWs. Here, the A valleys are 
degenerate, and for most barrier compositions the us- 
able energy range is limited by the A valley offset at the 
Si/SiGe interface. The system is less sensitive to strain, 
and the maximum usable energy range of 185 meV is ob- 
tained when the barriers have a Ge fraction of around 
89%. However, our designs in section ITVl use lower barri- 
ers with a Sio.4Geo.6 alloy composition, in order to obtain 
realistically wide layer widths, as described above. This 
composition provides a usable energy range of 90 meV. 

Fig. [ ^c)| shows the minima for (OOl)-oriented 
GeSi/Ge/GeSi QWs. In structures with barrier Ge frac- 
tions greater than around 75%, the L valleys form the 
bottom of the QW in the Ge layer. However, for lower 
Ge alloys, the A2 valleys are lowest in energy. The max- 
imum usable energy range of approximately 90 meV is 
obtained when the barriers have a Ge fraction of around 
0.85. It is worth noting that A2 states in the thin GeSi 
layers of QCLs will have confinement energies well above 
the A2 band edge. It may, therefore, be possible to ob- 
tain a larger usable energy range by using a lower barrier 
Ge fraction. 



C. Self-consistent Poisson-Schrodinger solution 

Self-consistent solutions of the one-dimensional time- 
independent Schrodinger equation and the Poisson equa- 
tion were found for the structures considered in this work, 
using a similar approach to those in Refs. [33| and [34]. 

The charge density over the length of the structure 
takes the form p(z) — e[d(z) — p e (z)], where e is the 
electron charge, d(z) is the ionized donor profile, p e (z) is 
the electron density profile and z is the spatial position. 
It was assumed that all donors are ionized, and hence d(z) 
is equal to the dopant distribution. An initial estimate of 
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FIG. 2: (Color online) Conduction band minima in the well and barrier regions of Si-based QWs relative to the 
bottom of the well. Energies are plotted as a function of the barrier alloy composition, and results are shown for 
structures grown on various Sii— x ,Gea:, virtual substrates. The figures show band minima in QW heterostructures 
comprising[(a71(001)-SiGe/Si/SiGe,[(b)1(lll)-SiGe/Si/SiGe, and \(cj\ (001)-GeSi/Ge/GeSi. 



the space-charge effect on the Hamiltonian was generated 
by solving the Poisson equation for a uniform electron 
distribution, in which p e {z) ~ N2d/L p , where N2D is the 
sheet doping density across a structural period of length 
Lp. 

A one-band parabolic effective mass approximation 
was used for the Schrodinger equation, which is justi- 
fied by the small confinement energy, and the large en- 
ergy separation from other energy bands. The quantiza- 
tion and density-of-states effective masses (m q and m-d 
respectively) were calculated for each material and crys- 
tal orientation using the method in Ref. [32| and are shown 
in table HI 

Intervalley mixing splits the A2 subbands in (001) 
Si/SiGe and the L subbands in Ge/GeSi heterostruc- 
tures. It is, however, only possible to include this ef- 
fect in complex heterostructures via computationally- 
expensive atomistic approaches such as tight-binding 35 
or pseudopotential calculations^ We have previously 
shown that the splitting energies are small in structures 
wider than 2-3 nm^ and have, therefore, omitted the 
effect in the present work. 

Three periods of the conduction band potential profile 
were used in the Hamiltonian for the QCL simulations, 
with box boundary conditions. This introduced an un- 
realistic limit to the spatial extent of the wavefunctions 
which were localized near to the boundaries. To elim- 
inate this effect, we replaced the Schrodinger equation 
solutions whose wavefunctions lay in the left and right 
periods of the structure with translated copies of the cen- 
tral period solutions. 

Having solved the Schrodinger equation, the subband 
populations, rij, were calculated as described in the next 
section, and a new estimate of the charge distribution 
was generated, using 

p e (z) = n val ^n,-|i/>j(z)| 2 , (11) 

i 

where n va i is the valley degeneracy and ipj is the wave- 
function of state j. The Poisson and Schrodinger 



equations were then solved iteratively to find the self- 
consistent solutions. 



III. CARRIER TRANSPORT MODEL 

We used a computationally-efficient scattering calcula- 
tion and a semi-classical approach to electron transport 
in QCLs as described in our previous worki 1 Simi- 
lar Boltzmann or rate equation based models have pre- 
viously yielded good agreement with experimental data 
for III-V mid-infrared 3 ^ and THz QCLsi Indeed, our 
carrier-transport model has been shown to calculate the 
current density and maximum operating temperature 
quite accurately for a 4.4 THz, 7- well chirped-superlattice 
GaAs/AlGaAs QCL that is similar in structure to the de- 
vices considered in section II VI 39 Although recent work 
on III-V QCLs has focused on coherent transport ef- 
fects, using nonequilibrium Green's functions 4 ^ or den- 
sity matrij* 4 ^— approaches, the semi-classical approach 
combines the computational speed and flexibility re- 
quired for extensive exploration of the parameter space of 
possible device specifications. We have assumed that co- 
herent transport effects are less significant in the bound- 
to-continuum devices considered in this work than in 
resonant-phonon QCLs, owing to the reduced thickness 
of the injection barriers through which electrons tunnel 
into the active region. Furthermore, the absence of reso- 
nant LO-phonon scattering may lead to longer dephasing 
times for coherent transport in group-IV materials than 
in III-V materials. 

As in our previous work ) 17 ' 37 our model includes elastic 
intravalley scattering due to interface roughness (allow- 
ing arbitrary interface geometries),^ alloy disorder ; 47 ' 48 
ionized impurities^ electron-electron interactions 5 ^ and 
deformation potential scattering for electron-acoustic 
phonon interactions? 5 ^ Intravalley optical phonon inter- 
actions are forbidden in A valleys due to the symmetry 
of the system^ but were included for L valleys via a 
zero-order deformation potential model. 51 ' 53 
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FIG. 3: Potential profile and electron probability 
densities for optimized (001) Ge/GeSi QCL designs, 
with emission frequencies and operating biases of |(a)| 
3.0 THz, 3.3kV/cm and [(b)] 3.8 THz, 3.5kV/cm. The 
upper and lower subbands involved in the radiative 
transition are shown as solid-bold and dashed-bold lines 
respectively. 



Intervalley phonon scattering was also described us- 
ing the zero-order deformation potential model, with the 
rates multiplied by the number of equivalent destination 
valleys. In L — > L scattering, all three destination val- 
leys are degenerate, and separated by a wavevector of 
the same magnitude. A phenomenological approach, de- 
scribed in Ref. [54| was used to describe the L — » L in- 
teractions by treating the combined scattering from all 
phonon branches as a single interaction. A — » A scatter- 
ing interactions are categorized as either g type, in which 
the destination valley lies on the same crystallographic 
axes as the source, or / type in which the destination 
valleys lie on a different axis. Phonon energies and defor- 
mation potentials for g and / interactions with longitudi- 
nal/transverse optical (LO/TO) and acoustic (LA/TA) 
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FIG. 4: Gain spectra for devices emitting near (a) 
3 THz and [(b)] 4 THz. 



phonon branches were taken from Ref. l55l The high- 
energy g-hO, /-LA and /-TO phonon interactions were 
determined using the zero-order deformation potential 
model. The lower-energy g-TA, g-LA, and f-TA inter- 
actions have no zero-order component in their deforma- 
tion potential, owing to symmetry selection rules, and 
were instead determined using a first-order modeL 56 i 57 
The bandstructure calculations in section [H] show that 
A4 quantum wells correspond to A2 barriers in (001) 
Si/SiGe heterostructures. The small spatial overlap of 
wavefunctions leads to very small A2 — ► A4 scattering 
matrix elements, and / transitions were therefore omit- 
ted in our model of (001) devices as a first-approximation. 

The steady-state populations rii for each subband were 
found using a rate-equation approach^ Intrasubband 
scattering rates were typically calculated to be an or- 
der of magnitude faster than intersubband scattering. It 
was therefore assumed that electrons settle between in- 
tersubband scattering events to a quasi-thermal Fermi- 
Dirac distribution. The distribution for each subband 
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was described by a quasi-Fermi energy Ep.i and a global 
electron temperature T e and the total subband popula- 
tion was found using 



80 



m = p 



2D k B T e \n 



1 



E F,i( T e) 



where p 2D = md/(iTfi ) is the two-dimensional density- 
of-states. 

A root-finding approach was used to determine the 
steady-state electron temperature at which no net gain 
or loss of kinetic energy occurred within the QCL, using 
the expression^ 



dE k 
dt 



EE' 



i(T e )E if Wif(T e ) = 0. 



(13) 



Here, Wit is the average intersubband scattering rate be- 
tween a pair of subbands i and j , summed over all scat- 
tering processes, and Eif is the energy difference between 
the subband minima. In the case of inelastic processes, 
the transition energy was modified as Eif —>Eif± huj q to 
account for the absorption or emission of a phonon with 
energy Uw q . 

The current density was estimated by considering the 
average scattering rates and the change in electron posi- 
tion for all intersubband transitions, 17 



J =77E^<a 1 E(^)/-^)^/' ( M ) 
P * f 



where L p is the length of a structural period of the QCL, 
n va i is the number of equivalent initial valleys and (z) is 
the expectation position for an electron in a given sub- 
band. 

The optical gain per unit length was calculated using 
G(lj) — cr(o;)/(£ocn r .) ) 38 ' 58 where n r is the real part of 
the refractive index of the active region stack and a is 
the real part of the optical conductivity. This is given by 
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(15) 

where is a lineshape function and sgn(-) represents 

the sign-function. The oscillator strength is given by 



Si 



2(m q m 2 )z 2 



(16) 



where, = is the dipole matrix element. A 

Lorentzian lineshape was assumed, with a linewidth of 
2meV, as is typical for the lasing transition in GaAs 
based THz QCLs^ 



IV. DEVICE DESIGNS 

Equivalent 3 and 4 THz QCLs were designed for each 
of the three material configurations, in order to simulate 
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FIG. 5: (Color online) Peak gain as a function of lattice 
temperature for devices emitting near 3 THz (red lines) 
and 4 THz (black lines) . The blue line represents an 
indicative figure of 31 cm^ 1 for the threshold gain. 



the relative performance of each system. To obtain a 
fair comparison between the materials, all designs were 
generated using an automated process. In principle, it 
would have been possible to generate slightly better de- 
signs than those presented in this section by manually 
adjusting the device structures at the start and end of 
the automated design process. However, we chose not to 
apply any manual design optimization in this work, in 
order to ensure that a consistent, reproducible and sys- 
tematic exploration of the design parameter space was 
used for each material configuration. 

Pure Si or Ge was used as the well material in order 
to prevent depopulation of the upper laser level via alloy 
disorder scattering. For each material configuration, a 
SiGe alloy was chosen for the barriers to maximize the 
usable energy range (as described in section HI)) , without 
introducing mechanical instability. 

A recent seven- well bound-to-continuum (BTC) (001) 
Ge/GeSi design was selected as a template for all the de- 
signs considered. This device gives a large simulated gain 
(~ 50 cm -1 ) at 3.5 THz, with a threshold current density 
around 300 A/cm 2 , an operating bias of 3.6 kV/cm^i and 
a maximum operating temperature of 136 Ki^ In this 
structure, doping was spread evenly over four wells and 
three barriers in the injector region of each QCL period, 
with a total sheet density of 8 x 10 10 cm~ 2 . Complete 
donor ionization was assumed. A similar seven- well BTC 
device has been demonstrated in the GaAs/AlGaAs ma- 
terial system, with an emission frequency of 3.66 THz at 
an operating bias of 4.15kV/cm and with a threshold 
current density of ^200 A/cm at low operating temper- 
atures. This GaAs/AlGaAs device was shown to have a 
maximum operating temperature of 116 K — the highest 
reported for a BTC THz QCL£ 

An automated design algorithm 6 ^ was used to generate 
a pair of new QCL designs from the Ge/GeSi QCL tern- 
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FIG. 6: (Color online) Relationship between electron 
temperature and lattice temperature for devices 
emitting near 3 THz (red lines) and 4 THz (black lines) . 



plate, with emission frequencies of 3 and 4 THz. In this 
method, the gain spectrum was calculated using a range 
of external electric fields from 3 to lOkV/cm. The thick- 
ness of each barrier and well was adjusted sequentially to 
maximize the gain in a frequency window of ±200 GHz 
around the desired emission frequency. The entire pro- 
cess was iterated until the algorithm converged on an 
optimal layer structure and field. The optimal virtual 
substrate composition was calculated for each device^ 
to ensure zero net stress across each period of the QCL. 

Equivalent 3 and 4 THz QCL designs were generated 
using (001) and (lll)-oriented Si/SiGe configurations. 
In each case, the thickness I of each layer in the tem- 
plate QCL structure was transformed according to I' — 

l^Jm q /m' q r& where the prime notation denotes param- 
eters of the Si/SiGe system in the appropriate orienta- 
tion. This transformation yields subband spacings ap- 
proximately equal to those in the original Ge/GeSi de- 
sign template. The automated design algorithm was then 
applied as described above. 

Parameters for each of the final QCL designs are sum- 
marized in table IIHI It can be seen that the Si/SiGe 
device designs generally require thinner layers than the 
Ge/GeSi designs, owing to the difference in effective 
mass. Epitaxial growth of QCLs in Si/SiGe may, there- 
fore, be more challenging. The total length of an ac- 
tive region period is also lower in Ge/GeSi designs than 
in Si/SiGe, which leads to a lower operating bias. The 
bandstructure and electron probability densities for the 
Ge/GeSi designs are plotted in Fig. [3J 



V. SIMULATED DEVICE PERFORMANCE 

The simulated gain spectra for all devices, operating at 
their respective design biases are shown in Fig. [4] It can 
be seen that the Ge/GeSi designs yield the highest gain 



at both frequencies, with peaks of 63 and 82 cm -1 for the 
3 and 4 THz designs respectively. The peak gains for the 
Si/SiGe designs are significantly lower: 25 and 38 cm" 1 
for the (lll)-oriented devices, and 5 and 3 cm -1 for the 
(OOl)-oriented devices at 3 and 4 THz respectively. We 
have previously calculated a threshold gain of 31 cm -1 
for a 15-/im-thick (001) Si/SiGe QCL active region in 
a copper double-metal waveguide structured By taking 
this threshold as an indicative figure, we predict that net 
gain is achievable for both Ge/GeSi devices, and for the 
4 THz (111) Si/SiGe device. 

The peak gains in the spectra decrease as the lattice 
temperature increases, as shown in Fig. [SJ owing to the 
reduction in population inversion. This is caused by 
a number of thermal effects, including electron leakage 
from the upper laser level via phonon emission, and by 
thermal backfilling of the lower laser level. Net gain is 
predicted for the Ge/GeSi devices up to lattice tempera- 
tures of 179 and 184 K for 3 and 4 THz emission respec- 
tively. The 4 THz (111) Si/SiGe device is predicted to 
yield net gain up to a lattice temperature of 127 K. 

The simulated temperature of the electron distribu- 
tion T e is plotted as a function of lattice temperature 
T in Fig. [5] At high lattice temperatures, T e is a lin- 
ear function of T and is approximately independent of 
bias. At low lattice temperatures, however, T e is deter- 
mined principally by the applied electric field. In the 
case of (001) Si/SiGe devices, the bias is relatively large 
(> 7kV/cm), and electrons therefore scatter preferen- 
tially into high-energy states. This yields high steady- 
state electron temperatures of 184 and 189 K for emis- 
sion at 3 and 4 THz respectively at a lattice temperature 
of 4K. The electric fields are lower in (111) Si/SiGe and 
(001) Ge/GeSi devices, owing to the greater lengths of 
the active regions. This leads to correspondingly lower 
electron temperatures of 127 and 129 K for (111) Si/SiGe 
devices, and 93 and 100 K for Ge/GeSi devices emitting 
at 3 and 4 THz respectively. The effect of thermal excita- 
tion upon device performance is illustrated in Fig. [7] It 
can be seen that the gain decreases monotonically as elec- 
tron temperature increases, owing to the thermal back- 
filling of the lower laser level. Ge/GeSi devices are able 
to operate with the lowest electron temperatures, and 
hence achieve the highest peak gains. 

The current density was calculated at the design 
bias for each of the devices. Current densities of 270 
and 380Acm~ 2 were predicted at the design bias for 
Ge/GeSi devices operating at 3 and 4 THz respectively. 
In Si/SiGe, current densities were calculated as 430 and 
460 A cm -2 for the (lll)-oriented devices and 210 and 
240 Acm~ 2 for the (OOl)-oriented devices at 3 and 4 THz 
respectively. The low operating currents in (001) Si/SiGe 
devices were due to the very low scattering rates, which 
result from the high A2 valley effective mass. The ra- 
tio of peak gain to current density was calculated as 
a figure of merit for each device at its design bias. 
Ge/GeSi devices were found to have the highest values 
(240 and 210cm/kA) followed by (111) Si/SiGe (57 and 
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TABLE III: QCL design parameters for each of the devices designed in this work, where fo is the emission frequency 
in THz, x s denotes the virtual substrate Ge fraction and F is the operating bias in kV/cm. Bold text in the layer 
structure represents barriers, while normal weighted text represents wells. Doped layers are underlined. 



Material 


fo 


Layers [nm] 


x 3 


F 


(001) Si/Sio. 8 Ge .2 
(111) Si/Si .4Ge .6 
(001) Ge/Geo.ssSio.is 


2.9 
3.1 
3.0 


2.1/0.8/6.4/0.7/3.5/1. 0/3.3/1. 2/3.1/1.4/3.0/1.4/2.9/1.6 
3.2/1.0/9.2/0.8/5.0/1.5/4.5/1.9/4.4/2.3/4.4/2.4/4.2/3.3 
6.7/1.2/15.1/1.4/11. 1/1.5/9.5 /1.8/8.7/2.3/7.7/3.5/7.1/4.3 


0.048 
0.146 
0.969 


7.1 
6.9 
3.3 


(001) Si/Sio. 8 Ge .2 
(111) Si/Sio. 4 Ge .6 
(001) Ge/Geo.gsSio.is 


4.1 
4.0 
3.8 


2.2/0.7/6.1/0.8/4.3/1. 0/3.1/1.0/3.1/1. 2/3.0/1.4/2.7/1. 5 
3.1/1.0/8.7/1.3/5. 3/1. 6/5. 0/2. 0/4.6/2. 1/4.0/2. 4/4.0/3.7 
5.8/1.0/15.3/1.4/12.3/1.6/9.9 /1.9/8.3/2.4/7.8/2.9/7.0/4.3 


0.045 
0.154 
0.970 


7.3 
6.9 
3.5 



(001) Ge/GeSi 
(111) Si/SiGe 
(001) Si/SiGe 




L00 150 200 250 

Electron temperature [K] 



300 



FIG. 7: (Color online) Relationship between peak gain 
and electron temperature for devices emitting near 
3 THz (red lines) and 4 THz (black lines). 



84cm/kA), and (001) Si/SiGe (25 and 14cm/kA) at 3 
and 4 THz respectively. We should note that our simula- 
tions of (001) Si/SiGe QCLs do not include A 2 ->■ A 4 in- 
tervalley scattering events, which would further degrade 
the predicted performance. However, as these structures 
already appear to be poor candidates for laser design, 
a more comprehensive transport model was considered 
unnecessary. Threshold current densities were calculated 
at T = 4K as 440, 210, and 330 A cm" 2 for the 4 THz 
(lll)-Si/SiGe, 3 THz Ge/GeSi and 4 THz Ge/GeSi de- 
vices respectively. 



VI. CONCLUSION 

We have presented a comparison between the sim- 
ulated performance of Si-based QCLs using the (001) 
Ge/GeSi, (111) Si/SiGe, and (001) Si/SiGe material con- 
figurations. A semi-automated design optimization al- 
gorithm was used, in order to provide a fair compari- 
son between equivalent designs. Our results show that 
(001) Ge/GeSi is the most promising system for devel- 
opment of a Si-based QCL. Firstly, the bandstructure 
calculations in section [H] show that the (001) Ge/GeSi 
and (111) Si/SiGe systems offer a >90meV energy range 
for QCL design, compared with only ^5meV in (001) 
Si/SiGe systems, owing to the large energy separation be- 
tween conduction band minima. This reduces the prob- 
ability of current-leakage via intervalley scattering, and 
allows a wider range of emission frequencies to be tar- 
geted. Secondly, the low L valley effective mass was 
found to yield a relatively long period length for the QCL 
active region. This reduces the operating electric field, 
and hence the current density and the temperature of the 
electron distribution. Net gain was predicted for both of 
the Ge/GeSi devices, but only one of the four optimized 
Si/SiGe devices. Ge/GeSi bound-to-continuum QCLs 
were predicted to operate up to temperatures of 179 and 
184K at 3 and 4 THz respectively, while the 4 THz (111) 
Si/SiGe device was predicted to operate up to 127K. 
These figures may potentially be improved via waveg- 
uide design optimization to minimize losses, or through 
the use of a resonant-phonon active region design. Never- 
theless, the predicted values exceed the highest-recorded 
operating temperature of 116 K for a 3.66 THz seven- well 
III-V BTC device^ 
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